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Abstract 

In these proceedings, we highlight recent developments from both theory and experiment related 
to the global description of matter produced in ultra-relativistic heavy-ion collisions as presented 
during the Quark Matter 2012 conference. 



1. Introduction 

Describing the matter produced in ultra-relativistic heavy-ion collisions with a limited set of 
global variables is a tantalising task. In fact, the challenge is not only a description of the matter 
created under such extreme conditions but also an understanding of the details of its evolution. 
It is fair to say that stunning progress has occurred in the last couple of years, as it is reflected 
in the latest results presented during this edition of the Quark Matter conference. A significant 
fraction of the discussions focused on the shape of the initial energy density of the collision in 
terms of fluctuations. Such a picture is now constrained by extremely precise measurements of 
the azimuthal anisotropy of emitted particles as discussed in Sec. [2] A fairly detailed "standard 
model" for the dynamical evolution of a heavy-ion collision emerges and is explained in Sec. [3] 
Recently, lattice QCD calculations have made tremendous progress in eliminating systematic 
uncertainties. As explained in Sec. [4] various groups now agree on the value of the critical 
temperature for the chiral transition within systematic and statistical uncertainties, and on the 
value of the interaction measure at high temperature to within 25%. Section [5] is dedicated to 
recent experimental estimates of chemical and thermal freeze-out variables: statistical thermal 
analyses are performed on hadron abundances and comparisons are made for nuclei and hyper- 
nuclei as well. Baryon-to-meson ratios in the intermediate pj region are studied at RHIC and 
LHC in order to investigate hadronisation mechanisms involving parton recombination. Blast- 
wave fits to hadron px-spectra are used to extract thermal freeze-out conditions as a function 
of collision centrality and beam energies. We conclude our review of global observables and 
fluctuations with an outlook towards further key questions that need to be addressed in the near 
future. 

2. Constraints on fluid dynamics: azimuthal anisotropy 

The fluid-dynamical description for the bulk dynamics of heavy-ion collisions received a 
tremendous boost in interest with the measurement of elliptic flow at RHIC, showing that fluid 
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Figure 1: Harmonic flow coefficients for charged particles: left and middle panel as a function of pj for two different 
centrality classes 1 6], right panel as a function of centrality percentile | 5]. Calculations from Refs. [l|,|7||. 

dynamics could quantitatively describe the collective flow of matter and coining the paradigm of 
the Quark-Gluon Plasma (QGP) being the "most perfect liquid" ever created. Later, fluctuation 
measurements were perceived as a handle on both initial state as well as transport properties 
of matter produced (first as a strongly interacting QGP - sQGP - and then as a hadron gas). 
Currently, we are at the crossroads where state-of-the-art modelling of initial conditions meets 
extremely precise experimental measurements of fluctuations. This is spectacularly illustrated 
by the presentation of B. Schenke with the "real (conference) timeQ" matching of the ATLAS 
event-by-event (EbyE) flow fluctuations for the anisotropic flow coefficients v„ = 2,ia 10] by the 
Impact-Parameter Saturated (IP-Sat) Glasma modelling of the initial conditions |3J], followed 
by 3+1 -dimensional EbyE relativistic viscous fluid-dynamical evolution with MUSIC fll. One 
of the striking features is that both pj and centrality differential V„=2,3,4,5 distributions from AT- 
LAS H and ALICE |5J], respectively, are correctly reproduced with a single value of the shear 
viscosity over entropy density ratio rj/s = 0.2 (see Fig.Q]). Such a value is consistent with the 
quantification of the systematical uncertainties performed by M. Luzum and for which the "con- 
servative" range of 0.07 < rj/s < 0.43 is obtained [8]. Nevertheless, it is a major achievement 
that IP-Sat Glasma initial conditions manage to reconcile V2 and v„>3 contrarily to MC-Glauber 
or MC-KLN initial energy density distributions for LHC energies. Moreover, with ultra-central 
collisions (2%o) recorded and selected by CMS |9[] as well as with the studies of correlations be- 
tween measured event-plane angles, one can expect new challenges for models and other insights 
into the fluctuations of the initial energy density. The main conclusions from the extraction of v„ 
with the Beam Energy Scan (BES) data at RHIC are similar but bring further constraints: when 
comparing to models, a low rjl s is always favoured (i.e., close to 0.16) but the px-differential v'3 
measured by STAR seems to increase progressively up to ~ 30% with increasing beam energy 
11.5 -> 200 GeV JH- With the addition of particle identification (tt*, K ± , p, and p) 
up to /7j - 4 GeV/c at yfs?m — 200 GeV, PHENIX reports an ordering for v„ = 3 4 with the same 
pattern than what is seen for V2 and no beam energy dependence for px-differential v„=3,4 for each 
particle species in a 20-60% centrality interval 111 111 . Although the analyses require more statis- 
tics, a significant effort is made by the experimental collaborations at RHIC and LHC to perform 
V2 estimates with higher cumulants in order to isolate non-flow and fluctuation contributions. In 
parallel of the well-known fact that the px-integrated V2 increases by ~ 30% from top RHIC to 
LHC energies Jl^L IT3L fl~4ll . it is underlined that the px-differential V2 extracted with the 4-particle 
cumulant method stays close in the range 1 < < 3 GeV/c from lower BES RHIC to LHC 



The limited statistics at the time of the presentation was later improved, see Ref. Q- 
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Figure 2: v„/n^ 2 vs. KEj/tiq with n = 2-A (from left to right) for n*, K*, p+p, in central (0-50%) Au-Au collisions at 
V*iSF = 200 GeV Q. 



energies IU5I1 . 

Concerning the LHC studies, it is also reported that (i) for all collision centralities, the mea- 
surements of integrated vi and V3 do not change if one uses 4- or 6-particle cumulant methods 
(similarly to what was noticed for V2); (ii) although V3 shows also an increase with pj, and then 
saturates around pj - 3 GeV/c, its centrality dependence is smaller than the one of V2 (to be ex- 
pected if initial density fluctuations are indeed at the origin of V3) 111 611 . It is important to mention 
here the excellent agreement between different experiments when comparing v>2 obtained with 
the 4-particle cumulant method 11411 . 

Several possibilities were investigated for finding a universal scaling for produced hadrons as 
a function of their number of constituent quarks (n q ) : (i) as shown in Fig. [2] PHENIX presented 
Vn/n'q 2 vs. KEj/riq for n ± , K ± , p+p which seems to hold for several harmonics (i.e., n = 2-4 for 
Au-Au collisions at y«NN = 200 GeV) and several collision energies [11]; (ii) STAR reported 
that V2/n q vs. (mj - mo)/n q works for most identified particles and up to top RHIC energies 
except an intriguing 2-cr deviation for the <p meson at the lowest energies of yfs^m = 7.7 GeV 
and 1 1.5 GeV 111511 : (iii) no scaling seems to hold up to LHC energies [16]. 



3. State-of-the-art modelling of heavy-ion collision dynamics 

The "standard model" for the dynamical evolution of hot and dense strongly interacting mat- 
ter created in ultra-relativistic heavy-ion collisions is the following: (i) individual parton-parton 
collisions copiously create gluons and quarks. The initial points of production are computed 
within various approaches. Some possibilities are the standard Glauber model for nucleon- 

nucleon collisions, the KLN model fl7ll . or the IP-Glasma model HI [It] which combines the 

1 — I 

Impact-Parameter Saturation (IP-Sat) model for the initial nuclear wave functions 111811 with the 
classical Yang-Mills dynamics of the produced color fields ("Glasma") lfl9ll . In a Monte-Carlo 
implementation of these models, initial conditions fluctuate from event to event. For an illus- 
tration, we show the transverse energy density profile at a time t = 0.2 fm/c after the collision 
in Fig. [3] The IP-Sat model shows fluctuations on the smallest scales, followed by the KLN 
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Figure 3: Initial conditions in the Glauber, the KLN, and the IP-Sat model, from Refs. [ 1 



model. The Glauber model has the smoothest initial conditions, i.e., the smallest fluctuations. 
How these initial conditions evolve further and on which time scale thermalization is achieved 
is presently not clear. Calculations within 4 theory 12011 show that actually thermalization (i.e., 
establishing an equation of state) is achieved prior to isotropization (i.e., equality of transverse 
and longitudinal pressure), cf. left panel of Fig. [4] This may necessitate the use of anisotropic 
fluid dynamics 12111 . 

Commonly, one assumes a rather rapid approach to local thermodynamical equilibrium, so 
that fluid dynamics is applicable. The initial conditions for fluid dynamics are then specified on 
a hypersurface in space-time (commonly a constant proper time surface t - tq - const.). From 
then on, the evolution of the system is determined by the conservation laws for net-charge and 
energy-momentum, 

= , d^Ti™ = . (1) 

For a unique solution, causal and stable formulations of dissipative fluid dynamics require in 
addition dynamical equations for the dissipative components of jV and T^ y . For a systematic 
derivation of these equations from kinetic theory, see Ref. [22]. If the microscopic interaction 
rates drop below the macroscopic expansion rate of the fluid, a fluid cell will not be able to sus- 
tain (approximate) local thermodynamical equilibrium any longer; the cell "freezes out". This 



freeze-out process is commonly performed along a hypersurface of constant temperature 123 ] 



How to compute the single-inclusive particle spectra in the presence of viscous terms was also 



reported at this conference II24I1 . The subsequent evolution of the system is either performed 



via a microscopic "afterburner" which takes into account elastic (and possibly inelastic) colli 



sions between individual particles before they hit the detector (see, e.g. Refs. 125L 12611 ). or one 
simply assumes complete cessation of microscopic interactions. Then, one simply computes the 
single-inclusive particle spectra along the freeze-out hypersurface in order to compare with ex- 
perimental data. Despite the apparent crudity of the underlying assumptions, the latter procedure 
is remarkably successful in reproducing the measured harmonic flow coefficients for charged 
particles, cf. Fig. Q] where they are shown for a calculation within relativistic dissipative fluid 
dynamics (lluD with IP-Glasma initial conditions and a constant shear viscosity-to-entropy den- 
sity ratio of rj/s = 0.2. One observes nearly perfect agreement for all flow coefficients for all 
centrality classes, except the most central one where the elliptic flow coefficient somewhat 
exceeds the experimental values. 
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Figure 4: Left: evolution of various components of the energy-momentum tensor within <j> theory, from Ref. 1201 , When 
2Pt + Ph (blue) approaches the energy density (black), fhermalization is achieved, while isotropization happens when 
the transverse (green) and longitudinal (red) pressure become degenerate. Right: the interaction measure (e - 3p)/T 4 
from various 1QCD calculations [ 29] . 



4. New lattice QCD results 

As outlined in the former section, there remains a large uncertainty in the choice of initial 
conditions for fluid-dynamical calculations. Thus, there is little hope in being able to perform 
a reliable determination of the equation of state (and probably also of the transport coefficients) 
from experimental flow data. In this situation, lattice QCD (1QCD) remains the most important, 
and moreover quantitatively reliable, tool to provide microscopic input for fluid-dynamical cal- 
culations. There has been considerable progress in the past year: 1QCD calculations from the 
HotQCD and Wuppertal-Budapest collaborations now agree on the value of (pseudo-)critical 



temperature for chiral symmetry restoration: T c = 154 ± 8 (stat.) ± 1 (sys.) MeV 12711 vs. 



T c = 155 ± 3 (stat.) ± 3 (sys.) MeV 12811 . Regarding the equation of state, one of the most 
important (systematic) uncertainties that remain to be resolved resides in the interaction mea- 
sure. As can be seen from Fig. [4] (right), while the values of the HotQCD collaboration agree 
fairly well with those of the Wuppertal-Budapest collaboration below T c , they are systematically 
larger by about 25% above T c . 

Another important input in dissipative fluid-dynamical calculations are the transport coeffi- 
cients. One would hope that, in the future, 1QCD calculations of these quantities reach a similar 
precision as those for the equation of state. Until then, one has to rely on more phenomenological 
approaches to extract values for 77/ s and the other coefficients, see e.g. lf3QL f^lfl . 

Correlations between, resp. fluctuations of, various quantities are an important tool to learn 
about critical behavior near phase transitions. This conference has seen a plethora of new data 
for these quantities. What is most striking is that 1QCD calculations are now in the position 
to make firm predictions that can be directly compared to experimental data. Figure [5] shows 
#2 _ (6N 3 Q )/(N Q ) (left panel) and R® 2 = (Nq)/(8N%) (right panel), where N Q is the net charge 
and 5Nq = Nq - (Nq) its deviation from its average IT32T1 . R^ is fairly insensitive to the value 
of the baryochemical potential /ig, so it can serve as a thermometer to determine the (chemical) 
freeze-out temperature. Once the temperature has been extracted, R® 2 , being fairly insensitive to 
temperature, can be used to extract the value of Hb.IT at (chemical) freeze-out. 
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Figure 5: fiE as function of T (left) and R® 2 as function of fis/T (right), from Ref. (32h . 

5. Further constraints on chemical and thermal freeze-out 

Once the QGP cools down and chemical freeze-out occurs, hadron abundances are very close 
to the ones recorded with the experimental setups. Statistical thermal model analyses of the 
hadron yields have been successfully describing the measurements on a large range of beam 
energies with few global variables only: the chemical freeze-out temperature T c h, the baryo- 
chemical potential p% and the volume V of the fireball. The STAR analyses of BES data at RHIC 
offer the possibility to investigate not only the beam-energy but also the centrality dependence 
of these parameters I33L For the lowest energies of ysNN = H-5 and 7.7 GeV, the results 
of the fits for peripheral collisions seem to systematically depart from the phenomenological 
description of 1 GeV per hadron ll35ll . The parameter extrapolation from RHIC to LHC energies 
is straightforward, since one expects T C h to get closer to the deconfinement temperature T c of 
1QCD and fig to vanish with increasing beam energy. However, some tension appears with the 
ALICE data analysis [36]: proton and anti -proton yields are low with respect to calculations 
for a temperature of T c h = 164 MeV (see Fig. [6] left panel) which may reflect the presence 
of annihilation during the hadronic phase. The importance of corrections for feed-down and 
secondaries from interactions with the detector material is underlined together with the benefit 
of vertex detectors for this purpose [37]. Hypertriton and anti-hypertriton are measured not only 
at RHIC but also at the LHC. Using the BES data as well as the additional statistics of the 2012 
run, the STAR Collaboration manages to obtain the excitation function of the ^H +^H px-spectra 
and refines the estimate of the lifetime r with measurements in a larger decay length interval. A 
combined fit of 2010+2012 results gives: r = 138 ±£. ps |38[]. Although ALICE extracts clear 
signals for ^H, ^H and ^He, no H-dibaryon peak is observed and therefore, upper limits on the 
yields of this hypothetical particle are derived j39tl . 

The observation of baryon/meson ratios at intermediate pj being significantly higher for top 
RHIC energy Au-Au collisions as compared to pp triggered the development of models recom- 
bining quarks from the QGP phase to produce hadrons. STAR tried to isolate the possible onset 
of parton recombination looking at the hyperon production using the BES program and in par- 
ticular the Q/0 ratio vs. pj: again, a difference seems to be seen for the lowest beam energies 
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2 Both grand canonical and strangeness canonical ensembles are used within the THERMUS 1 34] code, which lead to 
opposite variations of the extracted parameter. 
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Figure 6: Left: comparison of hadron integrated yields to pion ratios at mid-rapidity for RHIC and LHC central collisions 
with thermal model predictions [36]. Center: comparison between the p/x ratio as a function of pi for 0-5% central Pb- 
Pb collisions at y^NN = 2.76 TeV and several models |41]. Right: Evolution of 7\i n vs. (fij) extracted with simultaneous 
blast-wave fits to n ± , K*, p, and p p-p-spectra for different centrality intervals and beam energies at RHIC 1 33]. 



11.5 GeV [40]. As illustrated by the p/7r pj -ratio as a function of the Pb-Pb collision 
centrality in Fig. [6](center panel), the enhancement observed at RHIC still holds at LHC energies 
and the most central value (0-5%) is qualitatively compatible with several models using hadroni- 
sation mechanisms different from mere fragmentation 114 ill . A mapping of the kinetic freeze-out 
temperature 7\in vs. the mean transverse velocity (J3t) parameters is obtained performing blast- 
wave fits on n ± , K ± , p, and p pj-spectra IT331, l36ll . Both beam energy and centrality dependences 
are investigated with collisions from y/s^m = 7.7 GeV to 2.76 TeV. The evolution of the es- 
timated radial flow is smooth with (fij) increasing from ~ 0.1 to 0.65 and simultaneously 7\in 
decreasing from ~ 135 to below 100 MeV (see Fig. [6] right panel). For the most central colli- 
sions at the LHC, fluid-dynamical models are in general in good agreement with the identified 
/?x-spectra when a hadronic phase is included, in particular with antibaryon-baryon annihilation. 
The p T -spectra corresponding to 100 events simulated with MUSIC J3l+UrQMD [|42l|43l] during 
the time of the conference are very close to the measurements ll36Tl . 



6. Future prospects 

In the very near future, we perceive it to be mandatory to scrutinize the new measurements 
and check the consistency of the experimental results. Models have to be validated in the energy 
range from -\fs = 7.7 GeV up to 2.76 TeV. Then, systematic studies to further constrain the 
initial-state fluctuations are necessary. The hope is to get a better handle from detailed EbyE 
harmonic flow analyses, including event-plane angles. In particular, ultra-central events may 
serve as an interesting testing ground for models of the initial-state fluctuations. 

On the theory side, further developments point into the direction of including thermal fluc- 
tuations into the fluid-dynamical modelling, see Ref. 14411 . Of particular importance here seems 
to clarify how such fluctuations enter the dynamical equations of transient (second-order) fluid 
dynamics. Other interesting developments are the consistent description of the dynamics of the 
chiral order parameter within a fluid-dynamical framework [45]. Another interesting task for 
the future is the development of chiral anomalous fluid dynamics in order to assess the chiral 



magnetic effect [46] and the charge asymmetry observed by STAR 114711 and ALICE 14811 



Finally, explaining possible changes of hadron yields after chemical freeze-out (for instance, 
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due to baryon-antibaryon annihilation) and reconciling the statistical thermal model with data 
also seems to be an important task for future studies. 
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